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We examine core-level X-ray photoemission spectroscopy (XPS), X-ray absorption near-edge spectroscopy 
(XANES) and X-ray emission spectroscopy (XES) in the Falicov-Kimball model by using the exact solution from dy- 
namical mean-field theory. XPS measures the core-hole propagator, XANES measures the absorption of X-rays 
when the core electron is excited to an unoccupied electronic state of the solid, and not emitted, and XES mea- 
sures the spectra of light emitted as electrons fill the core-hole state created via some form of X-ray excitation. 
These three spectra are closely related to one another and display orthogonality catastrophe behavior at T = 0. 
We illustrate how to evaluate these spectra efficiently at finite temperature, with a primary focus on the details 
of XANES. 
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1. Introduction 

X-ray-based spectroscopy is one of the most important tools used to understand condensed matter 
physics. Experiments in X-ray spectroscopy are quite sophisticated and can measure a large range of 
different properties. The basic idea, however, is quite simple. X-rays are shone onto a sample and one 
measures either photons or electrons that are emitted after the scattering. Depending upon what is mea- 
sured, in what energy range, and how the initial X-ray energy is tuned relative to different absorption 
edges, one finds a "zoo" of different probes that are used within X-ray scattering. Here we focus on three 
of these probes, assuming we are using monochromatic X-rays shining on our sample. The first is X-ray 
photoemission spectroscopy (XPS) of the core levels. In this experiment, the X-ray is tuned far enough 
above an absorption edge that the photoexcited core electron has enough energy that it rapidly escapes 
from the crystal (without interacting with it) and effectively leaves the sample as a free particle. One mea- 
sures the rate at which electrons hit a detector as a function of the final energy of the electron, integrating 
over all angles (if one measure this as a function of the angle, then one can determine angle-resolved pho- 
toemission spectroscopy, but this is usually not important for core-level spectra which have essentially 
no momentum dependence). In this sense, it is a photon in, electron out spectroscopy. The second is X-ray 
absorption near edge spectroscopy (XANES) which measures the attenuation of X-rays as they are shone 
onto a material as a function of the incident X-ray energy. The X-ray energy is tuned such that the pho- 
toexcited electron does not have enough energy to leave the crystal, and hence it is similar to an optical 
absorption measurement, but with a different frequency of light. To get a spectra, one needs to change 
the X-ray frequency and measure the absorption as a function of the color of the X-ray light, just like how 
one measures and optical conductivity (after processing reflectivity data to extract the absorption). If the 
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electron does leave the crystal but with low-enough energy that it interacts strongly with the surface, 
then one has X-ray absorption fine structure (EXAFS), which we will not discuss further here. The third 
is X-ray emission spectroscopy (XES), where one has a system with an excited core-hole, and we measure 
the spectra of the light emitted from the material as the conduction electrons fill the hole and cause the 
sample to fluoresce. It is the stimulated emission of the core hole back down to the ground state, which is 
enhanced due to the presence of the X-ray radiation, and hence this is also an elastic process where one 
measures light emitted at the same frequency as the incident light, but in different directions. 

One can contrast these probes with inelastic scattering like nonresonant inelastic light scattering 
(NIXS) and resonant inelastic light scattering (RIXS) which are coherent second-order processes simi- 
lar to Raman scattering at optical wavelengths. XPS is a single-particle process that can be analyzed via 
the sudden approximation, XANES and XES are first-order processes in the matter-light interaction and 
correspond to two-particle processes, while NIXS and RIXS are second-order processes in the matter-light 
interaction which correspond to four-particle correlation functions, and won't be discussed further here. 

Before we discuss the formal developments in detail, we need to define our model. We work with 
the Falicov-Kimball model 1 1 1, with an additional core-hole interaction. For simplicity, we work in a hole 
picture, where the fermionic operators are all hole rather than particle operators. The spinless version of 
the Falicov-Kimball model involves the interaction of itinerant holes (described by the creation operator 
d\ at site i) which can hop between nearest neighbors and localized holes (described by the creation 
operator fj at site i) which do not hop, but interact with the conduction holes via a Coulomb interaction 
U. When we add the core-hole (described by the creation operator h\ at site i), it interacts via an on-site 
Coulomb repulsion with the conduction holes and Qf with the localized holes. The Hamiltonian is 
then 

Je = -Z t ij d\d j -nZdJdi +E f Y J f}fi + ^ E d\dj]f i + h\\ 

ij i i i i 

+ Z(Qdd f i d i + Q f flf i )h]h i , (1.1) 

i 

i 

where tij is a real symmetric matrix that is nonzero only if i and j are nearest neighbors and it equals 
t* 12 \fd on a rf-dimensional hypercubic lattice (we will be taking the limit d — «■ oo and using t* as the 
energy unit). The symbol /i is the chemical potential for the conduction holes, while Ef and Ej, are the 
site energies for the localized hole and core hole, respectively. The term ^J*™ denotes the first term of 
the Hamiltonian which includes the kinetic energy of the conduction electrons, while J3f! oc includes all 
of the remaining terms at site i, which are the local terms. Since the localized hole number operators 
f! f. and h'.h, both commute with 36, they are conserved quantities, and hence this model maps onto 
a generalized Falicov-Kimball model with spinless conduction holes, and spin-one-half localized holes, 
which can be solved exactly within dynamical mean-field theory (DMFT) as described in detail in the 
review in reference (2). 

The retarded local Green's function for the conduction holes is defined via 

G r d (t- t') = -W{t- t') ({^(0,4^)}+)' (1 ' 3) 

where the hole operators are expressed in the Heisenberg picture d-(t) = exp(L3ft)d. exp(-L3f t) and 
the angle brackets imply a trace over all states weighted by the density matrix p = exp(-f)3f)/3 with 
3 - Trexp(— $3€) and P-l/T the inverse temperature. The symbols {•, ■}+ denote the anticommutator. 
Because we have time-translation invariance in equilibrium, the Green's function depends only on the 
time-difference, and hence it is convenient to take a Fourier transform and represent it as a function of 
frequency u>. In addition, we need to define the self-energy £(&)), effective medium GqUo), and dynamical 
mean field X{<a). This is done as follows: (i) the local Green's function can be expressed in terms of a 
Hilbert transform of the self-energy with a weighting function given by the noninteracting density of 
states p{e) = exp (-e 2 ) / \fn which is a gaussian for the infinite-dimensional limit of a hypercubic lattice 

G r . (<u) = f dep (e) \— ; (1 .4) 

J ld + p-H{(u)-e 
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(ii) the effective medium is then defined via the Dyson equation for the local Green's function 

GoM = t^t— (1.5) 
and (iii) the dynamical mean field is then 

X{<d) = u)+u — . (1.6) 

It turns out that the dynamical mean field plays a prominent role in determining the different X-ray 
spectroscopies. 

Note that we are always working in the limit where the number of core holes in the system is so 
small, that the density vanishes in the thermodynamic limit, hence, in determining the dynamical mean 
field, we can set all h\h, operators to zero and recover the spinless Falicov-Kimball model. The core-hole 
operators, however, enter in important ways when we calculate the X-ray spectroscopies. Details for how 
to numerically solve for the Green's functions, self-energies, and dynamical mean fields of the spinless 
Falicov-Kimball model can be found in reference 0. 

In the next section we discuss the many-body correlation functions for the different spectroscopies. 
In Section 3, we describe how to solve exactly for these correlation functions within the DMFT approach. 
In Section 4, we show some numerical results, and we conclude in Section 5. 



2. Many-body formalism for the different correlation functions 

We start with a discussion of XPS. This is a photoemission process. If we invoke the sudden approxi- 
mation and ignore the momentum dependence of matrix elements which connect the free electron states 
to scattering states, then the XPS signal is simply given by the greater core-hole Green's function (see, for 
example, reference |5 |). This Green's function can be found from the so-called contour-ordered Green's 
function which is defined for each time lying on the Kadanoff-Baym-Keldysh contour (3]|4) depicted in 
figure [l] The contour starts at time t — 0, runs out to time t — t, then runs backwards in time and finally 
moves along the imaginary axis a distance of /3 = II T. The lesser Green's function has the first time lying 
on the upper branch and the second time lying on the lower branch, so the second time is ahead of the 
first on the contour. 




■i 3 



i 



Figure 1. The Kadanoff-Baym-Keldysh contour. The contour starts at time 0, moves along the real time 
axis to time t, moves back along the real axis to time 0, then extends down the imaginary axis to imaginary 
time -i/3. 



Before we define the contour-ordered core-hole propagator, we need to extend the definition of the 
dynamical mean field from real frequencies to times on the contour. This is done by defining the following 
function 

oo 

X c (t,t') = -- \ d«rm[AM]e~ k " (t ~ i ' ) \f{(a)-B c (i,i')], (2.1) 

7T J 

— oo 

where f(co) = l/(l + exp[/3(o]) is the Fermi-Dirac distribution function and0 c (f, t') is the generalization of 
the Heaviside unit step function to the contour, which equals one if i is ahead of t' on the contour, equals 
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zero when t is behind V on the contour and equals 1/2 when t-V. Note that on the contour, two times 
which numerically have the same value are only considered equal if they also lie on the same branch of 
the contour, and the designation ahead or behind refers to this ordering along the contour as illustrated 
by the arrows in figure [l] 

The contour-ordered Green's function for the core hole is then defined to be G c h {t,t') = 
-i(3' c h i {t)h j .{t')), where 3~ c denotes time-ordering along the contour. Within DMFT, a local Green's func- 
tion on the lattice can also be calculated as an impurity Green's function that evolves within the dynami- 
cal mean field, hence we can also write the core-hole propagator as 



G c h {t, = Tr 

-^imp 



e -A^ l0C j- c exp^ -i / dr / dfrf + (OA c (r,f')rf(f') \ h{t)hht') 



(2.2) 



Here, J€ la<L is the local Hamiltonian evaluated at one lattice site and 2i mp is the impurity partition func- 
tion calculated for an impurity evolving with respect to the local Hamiltonian and the dynamical mean 
field evolution operator, given by the exponential function above. We briefly discuss in the next section 
strategies that can be used to calculate the core-hole propagator. 

The XANES signal arises when the X-ray is tuned just slightly above the absorption edge, so a core hole 
is created and the photoelectron goes to an unoccupied state in the conduction band. When this occurs, a 
photon is taken out of the X-ray beam, so we measure this absorption by counting the number of photons 
that pass through the sample and comparing that to the number that pass when the sample is not there. 
Viewed in this way, the XANES experiment is similar to an optical conductivity experiment, except there, 
it is an optical photon causing an excitation of an electron from an occupied to an unoccupied state |6|. 
Nevertheless, the similarity between these two processes tells us that we can find the XANES signal by 
evaluating the imaginary part of the Fourier transform of the appropriate current-current correlation 
function. Here is how we find the corresponding current operator for a XANES experiment. First, we 
assume that the core-hole and the conduction electron are created at the same site, and for the moment, 
we assume the core hole will only create itinerant electrons, not localized electrons, hence the transition 
is to one band only Second, the photons have momentum q, frequency w q = cq, and a spatial profile for 
the vector potential that describes the light field given by an amplitude A multiplied by a shape function 
/q(r). We define the dipole matrix element at lattice site i via 

Af q (i) = <0|rfj— /q(r)e iq ' r e q - Vht\0), (2.3) 

which is the matrix element between a core-hole state h\\0) and an electronic state in the conduction 
band d^O), both of which need to be represented by their corresponding wavefunctions in coordinate 
space. Here e q is the polarization of the light field. In a translationally invariant system, the matrix ele- 
ment is independent of the lattice site i, so the current operator becomes 

j q = e-kVMqdjhJ + e iw i f M* h t d]. (2.4) 

The response of the current is given by the current-current correlation function. Noting that for XANES, 
we have no core hole in the initial state, we find the XANES spectra becomes 



|2 

h 



2\M I r 

t^ ANES (w a ) = — / Ate™* 1 (d\t)h{t)d{Q)h\Q)) , (2.5) 



where the local thermodynamic expectation value can be evaluated for the impurity problem, just like 
the core-hole propagator (implying we need to introduce the dynamical mean field evolution operator, 
etc.). 

The XES spectra is evaluated in a similar way, but this time there is a core hole present in the initial 
state, so different terms from the current-current correlation function contribute and we find 



2IM 1^ r 

x XES ( Wq ) = / Ate-™* 1 (d{t)hht)dhO)hm) ■ (2.6) 



o 
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In the remainder of this work, we will assume the matrix element M q is independent of momentum 
(and hence independent of the frequency <y q ). 



3. Formalism to evaluate the XANES correlation function 



We will focus primarily on showing the details for how to evaluate the XANES correlation function 
exactly. The XPS function can be evaluated using techniques similar to those developed earlier for the 
/-electron spectral function, which employ either a matrix-based formalism or a Wiener-Hopf sum- 
equation approach |7-9|. We are preparing a separate publication on this problem |10| which also ex- 
amines the numerical renormalization group approach 1 11 1 for the core-hole propagator. 

To calculate the XANES signal, we need to evaluate the four-point function that is defined by the 
functional derivative of the core-hole Green's function with respect to the dynamical mean field 



SX c {t 2 , h) 



sr c h{t)h\t')d(t l )dht 2 ]) + GUt, t')G c Ah, t 2 ) 



(3.1) 



where the symbol c denotes contour-ordered quantities and the itinerant hole Green's function G°, is 



defined analogous to Eq. 1 2.2 1 but with h[t)li*{t') replaced by d{t)d) It'); this Green's function can also 
be easily constructed from the retarded Green's function and the Fermi-Dirac distribution as a function 
of frequency because it is in equilibrium, but we do not give the details here. We need to evaluate the 
four-point function for t' = 0, t\ = and t 2 = t before we take the Fourier transform. 

Calculating this four-point function is straightforward, but tedious. To do so, requires us to define a 
fair amount of symbols before we can get to the final formula. In doing so, we follow closely the work in 
reference |9|, which shows how to represent these objects in terms of determinants of Toeplitz matrices, 
which can either be evaluated directly numerically, or can be solved in the large time limit asymptotically 
using Szegb's theorem I112H15I within the Wiener-Hopf approach (16] [17) as developed by McCoy and 
Wu (T$|. We don't provide details of those approaches here. 

We start by defining the fermionic Matsubara frequencies via ia> n = mT{2n + 1). Then, the impurity 
partition function can be represented as 2mq, = 2qq + Zqi + -Zio + Zn, with 

A(iw m ) 



j£ nhnj = e -P ( - E f n f +E h n h+Qfn f n h ) |j + e P()i-Un f -Q d n h )^ -Q 



1- 



, + \i-Un f -q d n d 



(3.2) 



with rif - 0, 1 and — 0, 1 the two possible fillings for each localized particle. Here \{iw m ) is the dynami- 
cal mean field evaluated at the Matsubara frequencies [which can be found by using the spectral formula 
given \ {(i>)}. Alternatively, these different partition functions can be defined over the contour via 



■2n;,nf — Tr 



e-^'^/^expi -ijdrj Aft 

I c c 



dht)X c [t, i')d(t') 



(3.3) 



so that the functional derivative satisfies 

S2n h rif 



Sk c {t2,t\) 



iTr 



pH i0 Hn h ,n f )^ exp \ -ij dij d? {i)\ c {t,t')d{?) ld(ti)rf t 



- --E nhri rGQd,n h nM\, h), 



(3.4) 



• l nf K J0d,ni 1 rif'. 

where God,n h n f is the effective medium projected onto the sectors of fixed and iif, and is defined in 
the equation above. 

Using these definitions, one has, for example, (f^f) = C2Toi + 2{ m -p. Since fiE^ — «■ oo for any 
reasonable temperature, we usually have 2\ nf — * (T is typically less than 0.1 eV, while Ef, is typically 
larger than 100 eV). 

Now we are ready to start the evaluation of the functional derivative. Note that we always have 
t > t' — 0, so we are evaluating the derivative of the greater core-hole Green's function, which has = 
in the initial state and the greater core-hole Green's function satisfies 



G>(t,0) = G> 00 (t,0) + G> m (f ) 0) 



Ml v 



(3.S) 
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with 



Gi m (t,0) = -i^-e- iE » t+c °°«\ (3.6) 

-^imp 

G Mi( f -°) = -i^e- i(£ " +( ?/ )t+c ° l(t) , (3.7) 
•2imp 

and 

C„ h „ f (t) =lndet c [f -<Q d {t,0)G odtnhnf ] . (3.8) 

The determinant is taken over the contour of the continuous matrix operator defined inside the brackets. 
Here, we have D c is the identity matrix on the contour, Q^(r,0) is a diagonal matrix equal to Q d on the 
top branch of the contour only [Q^(r,0) = QrfD [0,fl , with D |0,t| equal to the identity matrix for times on 
the upper branch, and equal to zero otherwise], and God,n h n^ is the continuous matrix operator for the 
effective medium projected onto the sectors of fixed nj, and rif. As it was already mentioned above, the 
greater core-hole Green's function {3.5} defines the XPS signal. 

The functional derivative of 2 nh n f is in Eq. 1 3.4 1. We also need the functional derivative of C nhTlf {t), 
which satisfies 

8X c {t z , fi) 

The functional derivative of the projected effective medium is straightforward to calculate and it yields 



f r r r 1 _1 r 8G 0dnn 1 

-^A[f-Q c d (t,o)G odinh „ f \ Q^ f -°) gAc(t2ttl) [- ^ 



8GQd,n h n f (t, r) _ 

~ Gad,n h n f it, tz)GQd,n h n f (h, t )■ (3.10) 



<5A c (r 2 ,fi) 

Substituting Eq. p35) into Eq. j3]9} yields 

8C„ h „Jf) I r ,-i 

g Ac(t2|fl) = -{ G od,B hW/ [» c -Qd(f,0)God, BfcB/ ] QSCf.OJGod,^^ ^ 

= G odtnftnf (.ti, t z ) - G^ dnhnf {ti,tz), (3.11) 

where we have defined 

G ?*« h n/ f i' t2) = | G 0d,n An/ [l C -Q^(f,0)G 0rf ,„ ;t „ / ] _1 }_ , (3.12) 
which satisfies the equation 



ht 2 



G 0d,n h n f ^> **> = G 0d , nh n f ih, fe) + Qd / dfGo d , Bfcn/ (fi, OG^^^ft, f 2 ). (3.13) 



Plugging this result into Eq. \3.l\ and taking the limits t\ = t' = and r 2 = t gives our final result 

X K ^(t) = (^ c Mt)dHt)dlO)hHo)) = -iGlw(tW (3.14) 



What is remarkable about this formula is that once one has projected onto the different static particle 
number sectors, the expectation values can be evaluated by factorization similar to Wick's theorem, even 
though we make no such assumption and the above result is exact. Since the core-hole propagator dis- 
plays the phenomena of the orthogonality catastrophe, where, in the limit as T — «■ 0, the spectral function 
diverges with a power law that depends on the coupling strength |6 19 1, we expect the same behavior to 
occur in the XANES (and XES) signals because they involve matrix determinants that are very similar to 
the XPS signal and should have the similar asymptotic behavior. 

In order to calculate this expression, we solve Eq. J3.13} via Cramer's rule, and then manipulate the 
determinants to get a uniform expression. We start by discretizing the interval [0, t] with N time slices 
(fi = and fjv - t), with a step size of Ar = tl[N - 1). Then we find we can approximate the XANES 
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correlation function by taking determinants of discrete (JV- 1) x fJV- 1) matrices (which can be scaled to 



(3.15) 

l-ArQ d G od>0 o(0) 



AfQrfGorf.oit^i - fjv-i) 
AfQdGod.oitfe ~ fjv-i) 

l-AfQ d G od) oi(0) 

Note that since each time t is independent of each other, this calculation can be easily parallelized; the 
t dependence in the determinant comes from which is nonzero only on the interval [0, t] and deter- 
mines the size of the matrix. The calculation for the XES signal is similar, but the details will not be given 
here. 



4. Numerical results 

For the numerical work, we will examine the Falicov-Kimball model on a hypercubic lattice at half- 
filling for the conduction holes and for the localized holes. We evaluate the XPS and XANES spectra for 
three different cases U = 0.5 which is a diffusive metal, U = 1 which is an anomalous metal that has a 
dip in the density of states at the chemical potential, and U - 2 which is in the Mott-insulating phase 
(due to space constraints we do not include XES results). We choose the core-hole interactions to satisfy 
Q d = Qf = 2. Both the XPS in Eq. |3^5j and XANES in Eq. |3.14| responses contain two contributions 
(groups of peaks), one from the sites with empty /-states and another from the occupied ones, separated 
by the Qf energy value and are represented by the blue dashed and red dot-dashed lines, respectively, in 
the figures below. 

For the cases of the diffusive (Figure [2} and anomalous (Figure [3| metals each of the contributions 
to the XPS response contain two peaks separated by the Q^ energy values. At high temperatures for the 
case of the diffusive metal (Figure |2^), we have in total four peaks (with two peaks merged due to the 
parameter selection = Qf). When temperature becomes lower (Figure^), the peaks become more 
pronounced with an almost Lorentzian profile and one of them becomes very sharp and at zero tem- 
perature transforms into the power law singularity of the Anderson orthogonality catastrophe |19|. For 
the case of the anomalous metal (Figure [3), the XPS response is almost the same with the appearance of 
additional shoulders at the peaks and now the edge singularity is much more narrow and sharper due to 
the decrease of the density of states at the chemical potential level, but it is still a power law divergence 
at zero temperature. The XANES response for these cases of the diffusive and anomalous metal contains 
two peaks, one of which transforms into a power law singularity at zero temperature. In other aspects, 
the behavior is similar to the XPS response. 

The behavior of the XPS and XANES spectra in the Mott-insulating case is completely different (Fig- 
ure [4f. The XPS response at high temperatures is similar to the one in the metallic phase, but when the 
temperature becomes lower, the first group of the peaks displays a gapped spectra similar to what is seen 
in the /-particle density of states [9]. The other group of the peaks displays a very narrow Lorentzian 
feature which changes at zero temperature from a power law singularity to a 5-peak. The XANES re- 
sponse also displays a sharp edge singularity, but now it is accompanied by a near edge structure which 
corresponds to the gapped Mott-insulator density of states. 



Af — » by repeating with different A t values). The formula is 



„XANES 



(?) = 



00_ e -iE h t 



imp 



x det 



Grjd.ooUi _ fiv) -AfQrfGorf,oo( f i _ h) 
God.oo ( - f;v) 1 ~ A tQd G d,oo (0) 



God.oo(Af) 

_ jfa -i(E h +Q f )t 



'imp 



x det 



Grjd,oi( f i _ f w) 
Grjd.oi^ ~ fjv) 



-AfQdGo(j,ooUiv-i ~ h) 



-ArQrfG rf,oi(fi - t z ) 
l-ArQ d G od ,oi(0) 



Grjdoi(A?) -AtQdGodfiiVN-i ~ h) 
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Figure 2. (Color online) XPS (a, c) and XANES (b, d) spectra for U = 0.5, Q d = Qf = 2, and T = 1 (a, b) and 
T = 0.3 (c, d). The black solid line gives the total spectra and the blue dashed and red dot-dashed lines 
correspond to the first and second contributions in the equations (3.5) and (3.14} , respectively. 





Figure 3. (Color online) The same as in Figure[2]but for U = 1. 
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Figure 4. (Color online) The same as in Figure[2]but for U = 2. 

5. Conclusions 

In this work, we have shown how to calculate three different X-ray spectroscopies in strongly corre- 
lated materials: XPS, XANES, and XES within the Falicov-Kimball model in the limit of large spatial dimen- 
sions. In all cases, one can find analytic expressions for these spectra in terms of determinants of con- 
tinuous matrix operators. These continuous operators can be discretized and evaluated on a computer 
using standard LAPACK routines. The XPS and XANES spectra have similar behavior, but the peak loca- 
tions and weights differ, indicating how they are measuring different processes. Generically, the XANES 
signal has fewer peaks, but the strong-coupling satellites are often appearing with higher weight. While 
these results do not incorporate complicated multiplet or spin-orbit coupling effects, which often are seen 
in real materials, they do include all of the strong-correlation effects in a simple model, allowing for an 
exact solution. 

A similar formalism to what was developed here can be used to solve for time-resolved XPS, XANES 
and XES experiments, where the system is pumped to a nonequilibrium state via a high-intensity fem- 
tosecond pulse, and then the system is probed with lower amplitude X-rays after some time delay with 
respect to the pump. The main difference in the formalism is that Green's functions lose their time- 
translational invariance and need to be solved using nonequilibrium DMFT |20|. We will present results 
for that work elsewhere. With the advent of newer high brilliance, short pulse width, X-ray sources (like 
the linac coherent light source, LCLS), we expect there to be more studies of pump/probe phenomena in 
strongly correlated materials, and it is possible probes like XPS, XANES, or XES, can shed light on inter- 
esting nonequilibrium dynamics. XPS has already been employed 1 21 1 to measure properties of the order 
parameter of charge-density-wave systems on ultrafast time scales, indicating that we are at the frontier 
of these types of experiments. 
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Tohhi/im po3B'5i30i< TeopiV flMHaiviiHHoro cepeflHboro nonn pj\5\ 

pi3H0MaHiTHMX f\OCJ\\f\\B 3 X-npOMeHflMI/1 B MOfle/li 

<t>a/i i ko Ba-Ki m 6a/ia 

A.M. LLlBaMKcP, flx.K. cbpipiKcP 

iHCTI/lTyT Cpi3M KM KOHfleHCOBaHHX CHCTeM HAH YKpaiHI/1, By/1. CBEHL4il^bKOrO 1, 79011 /lbBiB 

tt>i3WHHWM c))aKy/ibTeT, YHiBepcwTeT flxopflXTayHy, BaLuwmroH, OKpyr Ko/iyM6m 20057, CLUA 

flocniflxyiOTbai pemteHiBCbKi cneKTpw cfiOTOe/ieKTpoHHOi eiviiciT (XPS), nepeflKpafioBoro nor/WHaHHfi (XANES) 
Ta cfiOToeMiciT (XES) fl/w MOfle/ii ct>a/iiKOBa-KiM6a/ia BMKOpucTOByKDHi/i TOHHi po3B'H3KH Teopii fli/maMHHOro ce- 
peflHboro no/in. XPS BUMiproe nponaraiop flipKH mohhoto 3a/iniuKy, XANES — nor/inHaHHfl X-npoMeHiB, koiua 
e/ieKTpoH MOHHOro 3a/iniiiKy 36yflxyerbCH b He3anoBHeHy e/ieicrpoHHy 30Hy TBepfloro Ti/ia, i XES — cneicrp bh- 
npoMmeHoro CBiT/ia npn 3anoBHeHHi e/ieicrpoHOM Hanepefl CTBOpeHOiX-npoMeHflMH flipKH fiom-ioro 3a/iwLiJKy. 
Bci L(i Tpn ™nw cneKTpiB TiCHO n0B'fl3aHi Mix co6okd i npoaB/iaiOTb oco6/inBOcri Tuny KaTadpocfin opToroHa/ib- 
HOCTi npn T = 0. riOKa3aHO, hk ecfieKTHBHO po3paxoByBa™ TaKi cneKTpn a/w CKmneHnx TeMnepaTyp, 30KpeMa 
3 fleTa/ibHiLUWM Haro/iocoM Ha cneiapn XANES. 

lOuoHOBi c/iOBa: ipOToeneKTpoHHa peHTfeHocneKTpocKonin (XPS), peHTfeHOcneKTpocKonia nepeflKpadOBOro 
nor/MHaHHH (XANES), reopin fli/iHaMi^Horo cepeflHboro nom, MOfle/ib <S>aritKOBa-K\M6ana 
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